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Abstract. We have calculated first and second order corrections to several sum rules 
measured in deep inelastic lepton-hadron scattering. These corrections, which are 
due to heavy flavors only, are compared with the existing perturbation series which 
is computed for massless quarks up to third order in the strong coupling constant 
^1 as- A study of the perturbation series reveals that the large logarithms of the type 

^«0 . In* Q^/rn^ dominate the perturbation series at much larger values than those given 

by the usual matching conditions imposed on the as{n). Therefore these matching 
conditions cannot be used to extrapolate the running coupling constant from small /j, 
to very large scales like /x = Mz- An alternative description of the running coupling 

Q^ ■ constant in the MOM-scheme is proposed. 

^'- INTRODUCTION 

Oh! 

^ ■ Structure functions, measured in deep inelastic lepton-hadron scattering (see 

■ Fig. 1) 

X' hih) + H{p)^hih) + 'X', (1) 

H ; 

5r ! where 'X' denotes any inclusive final hadronic state, provide us with an excellent 

test of perturbative QCD. The in and outgoing leptons are represented by /i and 
I2 respectively and the hadron is denoted by H. On the Born level the reaction 
proceeds via the exchange of one of the vector bosons V of the standard model 
which are given by 7, Z and W^. The virtuality of the vector boson V and the 
CM. energy are given by g^ = —Q^ < and S = {p + kiY respectively. Further 
the scaling variables are defined by x = Q"^ /2p.q and y = p.q/p.ki. The structure 
functions show up in the polarized and unpolarized cross sections. Starting with 
the spin-averaged cross section for V = 'j we obtain 



^^ Talk presented at the International Workshop on Hadron Physics 'Effective Theories of Low 

Energy QCD', Coimbra, Portugal, 10-15 September, 1999. 

^^ Work supported by the EC network 'QCD and Particle Structure' under contract No. FMRX- 

CT98-0194. 




'X' 



FIGURE 1. Kinematics of deep inelastic lepton-hadron scattering. 
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For V = W^ and Q^ <^ M^ the unpolarized cross section for the charged current 
process becomes 
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Finally we are also interested in the spin structure functions which are measured 
in polarized scattering. In the case of l^ = 7 the cross section becomes 
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from which one can extract the longitudinal spin structure function gi. The trans- 
verse spin structure function (72 is obtained from the cross section 
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where is the angle between the spin of the hadron and the plane spanned by the 
lepton momenta ki. Further we have indicated between the brackets the spin of 
the incoming lepton /i and we have neglected in the expressions above all power 
corrections of the type M'^/S where M denotes the mass of the hadron H (see Fig. 

1). 

One of the most important predictions of QCD is the Q^-evolution of the struc- 
ture functions. However the a;-dependence cannot be predicted yet. The latter 



follows from non-perturbative QCD which is at such a premature stage that the 
x-dependence cannot be determined. This problem can be avoided when one inte- 
grates over the latter variable which leads to the sum rules discussed below. Notice 
that this integration requires a full knowledge about the range < a; < 1. Since 
there are no data available for x < 10"'^ and x > 0.8 one has to make extrapolations 
into this region. This will introduce an error which is very hard to estimate. In 
the past various techniques have been used to compute the sum rules. The most 
known among them are 

1. Infinite momentum frame techniques used in current algebra. 

2. Dispersion relations which are derived from the Regge behavior of the structure 
functions F{x, Q'^) at small x. 

3. Operator product expansion (OPE) techniques which in leading twist are 
equivalent to the parton model. 

Using the last technique the derivation of the sum rules proceeds as follows. First 
let us define the hadronic structure tensor which appears in the derivation of the 
cross sections presented above. It is given by 
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This tensor can be split up into a symmetric and a antisymmetric part i.e. 
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In the Bjorken limit {Q"^ -^ oo and x is fixed) the integrand in Eq. (6) is dominated 
by the lightcone z'^ = 0. This allows us to make a lightcone expansion of the current 
commutator 



iAz), -m] = y: e c^^'^(^v')o^-(/x^ o) , 



(9) 
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where for convenience we have dropped the Lorentz indices of the currents. In 
the expression above r and A^ denote the twist and spin of the operator O^'^ 
respectively. The latter and the singular coefficient function C^'"^ are understood 



to be renormalized where n represents the renormahzation scale. After insertion of 
expression (11) into Eq. (6) one obtains the A^th moment of the structure function 
which equals 

I' dxx^-^ Fix, Q^) = E (^) '" A(^)-(^^) C(^)- (^) . (10) 

In momentum space the operator matrix element (OME) and the coefficient func- 
tion are defined by 

^^^^•^(/^') = (p|O^'^(/^',0)|p), (11) 

and 

Cmr (91\ = |rf4^e''?-C^.-(^V) , (12) 

respectively. Limiting ourselves to twist two contributions only the non-singlet part 
(w.r.t flavor) of the structure functions is determined by the quark operator 

O^ = O^'"^-"" (x) = ^{x)Y'D^'^ ■ ■ ■ D^'^Xk^ix) , (13) 

where V^ denotes the covariant derivative and A^ are the generators of the flavor 
group SUpiuf). In the case of A^ = 1 the following operators are conserved 

0^(a;) = i){x)-i^\ki'{x) and , Of{x) = V5(a;)7VAfcV^(x) , (14) 

which means that they do not have to renormalized. Hence the corresponding 
OME's and coefficient functions are independent of the scale /i i.e. 

f'dxF{x,Q')=A^^^cl}\Q'), with (p|O,^(0)|p) = 4V, (15) 

Jo 

where A[. can be determined from SUp^nf) and a low energy theorem. In general 
the coefficient function C^-^^ can be expanded in a power series in Q;s(/i^). However 
in two cases this function becomes trivial. The ffist example is the Adler sum rule 

[1] 



/' - [Ff (x, Q') - Fr(a;, Q')] = K{nf) , 
Jo X ^ J 



A'(3) = 2 + 2 sin' ^e {SUf{3)) , K{4) = 2 {SUf{^)) ■ (16) 

Here 6c denotes the Cabibbo angle and for the constant K{nf) we have quoted the 
values given by the flavor group SUp{nf) for rif = 3, 4. The second example is the 
Burkhardt-Cottingham sum rule [2] given by 

1 
dxg2ix,Q^) =0. (17) 



The above sum rules above have the following properties 



1. The values on the righthand side are independent of the method used for the 
derivation of the sum rules. 

2. No power corrections of the type (A^/Q^)* (higher twist). 

3. No mass corrections (e.g. due to heavy flavors). 

4. No QCD corrections. 

The coefficient function has a nontrivial form in the following cases provided the 
results follow from OPE. The first one is the polarized Bjorken sum rule [3] 



^g,m^ / dx \gT{x,Q^)-gr{x,Q^) 
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The second one is represented by the unpolarized Bjorken sum rule [3] 
AFi(g2) = r dx \Frix,Q^)-F,'Pix,Q')] = K{nf)C^^{nf,Q'' 



K{3) = 1 + sin^ 9^ {SUf{3)) , K{4) = 1 {SUf{^)) , 
The third one is the Gross-Llewellyn Smith sum rule [5] 



(19) 
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F^^ix, Q') + F^ix, Q')] = K{nf) C^«(n^, Q' 



K{3) = 6-2 sin=^ 9^ {SUpiS)) , K{A) = 6 (5f/F(4)) . 
The coefficient function can be expanded in a^ as 

C%n,,Q')=± H''^^^'A \l{nf), r = g„F„F,. 
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The sum rules in Eqs. (18)-(20) have the following properties 

1. For rriq = the c[ are known up to order a^ for i < 3 (see [6], [7]). 

2. Higher twist corrections are estimated in [8] 

3. Heavy flavor corrections are computed up to order a^ in [9] 

In leading twist the coefficient functions can be computed in the QCD improved 
parton model. In this model the structure function can be written as 
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FIGURE 2. heavy flavor loop contribution to quark-quark scattering. 



where /^ denotes the bare parton density and the parton structure function Ta 
represents the QCD corrections including the Born approximation. Notice that 
Ta contains mass singularities. They are regulated by giving the light quark a 
mass which is sufficient in the non-singlet case. Subsequently these singularities 
are removed via mass factorization which reads as follows 

^a(^)=Aaa(^)®C,(^). (23) 

Hence the hadronic structure function can be written as 

i^(Q')=Ee^/f(/^')®C„(^), with /r(/i^) = ^aa(4)®/r' (24) 

where /^ denotes the renormalized parton density which like the coefficient function 
depends on the factorization scale [i. Performing the Mellin transformation one 
obtains 

'cia;a:^~iF(x,g2) = A(^)(^2)^(jv)(Q!)^ (25) 

[X^ 



so that one can make the following identifications 

H^ Jo n^ 

which yields the coefficient function C'-^^ 



A^''\f^') = fdzz^-' T.ellf{z,f^'), C(^H^) = f'dyy^-'C{y,%), (26) 
Jo „ u^ Jo u^ 



HEAVY FLAVOR THRESHOLDS IN THE STRONG 
COUPLING CONSTANT 

Coupling constant renormalization in the case of light flavors (massless quarks) 
conventionally proceeds in the MS-scheme. Using n-dimensional regularization the 



coupling constant gets a dimension i.e. g^ M*^~" where (7" is the bare (unrenor- 
mahzed) couphng constant. The renormahzation of the perturbation series is then 
carried out by the substitution 



fi-r -^ dsifJ'^ 



-^^/3o(n/) <^ - + 7i? - In 47r + In — 



(27) 



with Us = 9^/4:71 and f3o{nf) = 11 — 2/3 nj. The renormahzed coupling constant 
as satisfies the differential equation 

^' '^'^f'/^ =PMn,,^')), -> a.(n„/.^) = — i^ + .... (28) 
dfi (3o{nf)\nf2 

The solution obtained from the equation above is quite appropriate since the scale 
H is much larger than the light quark mass. However in the case of heavy quarks 
we encounter a problem which will be illustrated by the following example. Let us 
consider (light) quark-quark scattering as shown in Fig. 2. Including the heavy 
quark loop the cross section behaves like 

a ~ aliuf, ^2)i + 2al{nf, f^')]^^) + " " " , (29) 

where ni(x) denotes the one- loop vacuum polarization function which behaves 
asymptotically as 

This behavior leads to large corrections in the perturbation series which have to be 
resummed as follows. First we have to make the substitution 

asijif,!!^) -» as{nf + l,/i^)[l + ^^(n/ + l,/i2)/3o,Q In ^ ■ ■ ■] , (31) 

in Eq. (29) so that for —t ^ m^ the cross section behaves like 

1 I u^ 

a ~ al{n^ + l,/!^)- + 2al{nf + 1,;x2)-/3o,q In ^^ . (32) 

Next as{nf + 1,/i^) satisfies the massless renormahzation group equation in (28) 
but now for n/ + 1 light fiavors with the boundary condition [10] 

asiuf, Anf,n^) = asiuf + I, A„^+i, ^u^) , (33) 

where the matching value is given hj fi = m. In this way one obtains a better 
behavior of the perturbation series than shown by Eq. (29) in particular for n"^ = 
—t. However the representation of the cross section in (32) at n/ + 1 light fiavors 
is only correct if —t ^ m^ . This follows from the vacuum polarization function 
Hi(— 1/?72^) which gets close to the logarithm in Eq. (30) for —t > 21 m^ within 
10%. For —t ^ m^ one does not reach the asymptotic form of Hi(— t/m^) so that 
in this case it is much better to use the representation for the cross section at nj 
fiavors in Eq. (29). It also shows that the matching scale /i in Eq. (31) has to 
chosen at a much larger value e.g. /i^ = 21 rn^ instead of /i^ = 771"^. 



HEAVY FLAVOR CONTRIBUTIONS TO THE 
POLARIZED BJORKEN SUM RULE 

The sum rule in Eq. (18) receives contributions from the following subprocesses 
(see Fig. 2 in [9]) 

q + l*^q + Q + Q, 

q + '~f* —^q + g one-loop gluon self energy correction , 

q + 'J* -^ q two-loop vertex correction with one-loop gluon self energy . (34) 

The correction to the sum rule up to order a^ can be expressed as 



C^i(n;,g^) = l+ ^-"^^'"^ ai 
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where Hq denotes the heavy quark contribution. The two regions of interest are 

Q'«m'^HQ(^]^%ln%, (36) 

Q-'^^rr? ^ Hq (\\ = -a,Po,Q In ^ + 62 • (37) 

The first limit shows that heavy quarks of mass m decouple from the perturbation 
series when their mass is much larger than the virtuality of the photon. In the 
region Q^ ^ m^ we observe the same logarithmic behavior as for IIi in Eq. (30). 
However here C^^{nf, Q"^) becomes asymptotic for Q^ > 40 m^ within 15%. Hence 
the value of Q^ is about twice as large as the one found for —t in quark-quark 
scattering. This we have checked for charm production. Here we choose Uf = 3 
and m = irtf. = 1.5 GeV/c^ in Eq. (35) and the results are shown in table 1 (see 
[9]). Further we observe that the heavy flavor corrections to the sum rule (21) 



TABLE 1. QCD corrections to the polarized Bjorken sum rule. 




(52 (GeV/c)2 light flavors charm contribution^ charm contribution*^ 


charm contribution'^ 


3 0.795 -0.688.10-3 0.451.10-2 
10 0.883 -0.788.10-3 0.569.10-3 
100 0.931 -0.107.10-2 -0.0912.10-2 


-0.265.10-1 

-0.586.10-2 

-0.121.10-2 



^ Exact formula Eq. (35) . 

^ Asymptotic formula according to Eq. (37). 

'^Estimated order a^ correction according to [11]. 



are very small. They are even smaller than the estimated order a^ corrections in 
[11] in spite of the fact that the heavy flavor contribution Hq starts in order a^. 
The reason that the logarithmic behavior in process (36) emerges at a much larger 
scale than the reaction in Fig. 2 can be explained as follows. In Fig. 2 the gluon 
momentum k entering the self energy is an external variable since t = k^. However 
in the case of process (36) one has to integrate over k so that the latter is not an 
external variable anymore. Therefore the role of t is taken over by Q^. At very 
large Q^ the perturbation series in Eq. (35) can be improved by following the same 
procedure as outlined below Eq. (30). The result is that instead of Eq. (35) one 
obtains the improved perturbation series 
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which behaves much better for the choice /x^ = Q^ provided Q^ ^ rn^. However 
the expression above is not a correct representation when Q^ ^ rn^ as we have 



already seen below Eq. (33). Hence jj? 



m 



is not a good matching scale where 



the running coupling constant jumps from Uf to n^ + 1 flavors. Therefore this scale 
has to be chosen at a much larger value than what is usually done in the literature. 
This also means that extrapolations from a measurement of as{nf, /x^) at small fi to 
large values like /i = Mz have to be distrusted. A scheme (hereafter called MOM) 
which incorporates the effect of the heavy flavor thresholds in the running coupling 
constant in a much better way is proposed in [12]. Here one simply resums the 
vacuum polarization function as follows 
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where the coupling constant as{3,n'^) is represented in the MS-scheme. Expression 
(39) agrees rather well with the numerical solution of the renormalization group 
equation in Eq. (28) where the /9-function depends on m. This in particular holds 
when the two-loop self energy contributions H2 are included. From the results 
discussed above we have seen that the scales where both Hi and Hq behave log- 
arithmically are relatively much closer to each other than to fi = m. Hence it is 
better to put nj = 3 in Eq. (35) and to make the replacement 
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than to follow the procedure outlined below Eq. (30). In this way one gets a 
continuous transition through the flavor thresholds, also with respect to derivatives, 
of the perturbation series of C^^ in Eq. (35). Similar results were also found for the 
sum rules in Eqs. (19) and (20) (see [9]). Finally we conclude 

1. Heavy flavor corrections to the sum rules are very small and they do not exceed 
the estimated order a^ contributions due to light quarks. 

2. The matching scale in the literature (here fi = m) where the running coupling 
constant jumps from Uf to n/ + 1 flavors is chosen to be too small. 

3. The MOM scheme as chosen in Eq. (39) gives a better representation for the 
perturbation series than the usual MS-scheme with the matching condition in 
Eq. (33). 
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